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We investigated the low energy excitations in the parent compound NdFeAsO of the Fe-pnictide 
superconductor in the peV range by a back scattering neutron spectrometer. The energy scans on 
a powder NdFeAsO sample revealed inelastic peaks at E = 1.600 ±0.003^eV at T = 0.055 K on 
both energy gain and energy loss sides. The inelastic peaks move gradually towards lower energy 
with increasing temperature and finally merge with the elastic peak at about 6 K. We interpret the 
inelastic peaks to be due to the transition between hyperfine-split nuclear level of the ^*^Nd and 
^''^Nd isotopes with spin I = 7/2. The hyperfine field is produced by the ordering of the electronic 
magnetic moment of Nd at low temperature and thus the present investigation gives direct evidence 
of the ordering of the Nd magnetic sublattice of NdFeAsO at low temperature. 



PACS numbers: 75.25.+Z 

The recent discovery of superconductivity in doped 
iron-arsenide-oxideJ^ has heralded a new era in su- 
perconductivity researctP. After the first report that 
LaFeAsOi_a;Fa; became a superconductor below Tc — 
26 K, related compounds with higher critical temper- 
atures, up to 55 K in SmFeAsOi-xF^;, were quickly 
developecP. It is now accepted that these materials com- 
prise a second class of high- Tc superconductors, distinct 
from the superconducting cuprates discovered more than 
20 years agcP. The parent compound LaFeAsO crystal- 
hzes in the tetragonal ZrCuSiAs-type structure (Space 
group PA/nmm)^ . It is currently believed that the su- 
perconductivity in doped LaFeAsO is, like in cuprates, 
is not conventional BGS-type but intimately connected 
with magnetic fluctuations and with a spin density wave 
(SDW) anomaly within the FeAs layerd^. Undoped 
LaFeAsO undergoes a SDW-driven structural phase tran- 
sition around 150 K, in which the lattice symmetry drops 
from tetragonal to orthorhombic^ and anomalies in spe- 
cific heat, electrical resistance and magnetic susceptibil- 
ity occur. The antiferromagnetic ordering of the mag- 
netic moments (0.36/iB/Fe) was found below Tn = 134 
K by neutron scattering^. Electron doping by F or by 
O deficiency, as well as hole doping with Sr suppresses 
the phase transition and the tetragonal phase becomes 
s upercon ducting at in the temperature range 25-41 
K l^ l ^° l ^^ l Thus there is strong evidence that supercon- 
ductivity in LaFeAsO is mainly due to specific structural 
and electronic properties of the (FeAs)"^" layers. The mag- 
netic structures of the RFcAsO family have been investi- 
gated by neutron diffraction. As we mentioned before the 
Fe moments order below about 150 K accompanied by 
a tetragonal-to-orthorhombic phase transition at about 
the same temperature. The rare-earth moments order at 
lower temperatures. The magnetic ordering of NdFeAsO 



has been in vestig ated by specific heatl^^l and neutron 
diffractioiP^HHHI NdFeAsO undergoeJiS a tetragonal- 
to-orthorhombic structural phase transition at T5 « 142 
K. The Fe moments order antiferromagnetically in the 
stripe-like structure in the (a-b) plane below Tj\j ~ 137 
K but change from antiferromagnetic (AFM) to the fer- 
romagnetic (FM) arrangement along the c axis below 
« 15 K. The Nd moments order antiferromagnetically 
below about 6 K with the same structure. However neu- 
tron diffraction measurements suggest that the Nd mo- 
ments are polarised below T* 15 K. Here we report a 
high resolution inelastic neutron scattering investigation 
of the Nd magnetic ordering by studying the transition 
between hyperfine split nuclear levels of Nd nucleus gen- 
erated by the hyperfine field at the Nd nucleus due to the 
ordered electronic magnetic moment of Nd ions. This rel- 
atively less-known technique has been rece ntly d eveloped 
and tested in several Nd-based compound4^^*'25| 
elastic peaks observed in NdFeAsO below about Tjvd « 6 
K show directly the magnetic ordering of Nd ions at low 
temperature. From these studies we established that the 
energy of nuclear spin excitations is proportional to the 
ordered magnetic moment in Nd compounds. Here we 
used this linear relationship to estimate the ordered mag- 
netic moment of Nd ions in NdFeAsO. By this technique 
we only observe the inelastic signal due to the magnetic 
ordering of Nd ions because the spin-dependent neutron 
scattering cross section of ^^Fe is too small (0.5±0.1 barn 
for ^^Fe compared to 56 ± 3 barn for -'^"'^Nd ) to observe 
the signal due to the ordering of Fe moments. So the 
presence of inelastic signals in NdFeAsO at low temper- 
ature gives direct evidence for the magnetic ordering of 
Nd ions. 

The principle of the method of studying the hyper- 
fine interaction by the inelastic neutron scattering can be 
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summarizecp£12II as follows: If neutrons with spin s are 
scattered from nuclei with spins I, the probability that 
their spins will be flipped is 2/3. The nucleus at which 
the neutron is scattered with a spin-flip, changes its mag- 
netic quantum number M to Af ± 1 due to the conser- 
vation of the angular momentum. If the nuclear ground 
state is split up into different energy levels Ei\i due to 
the hyperfine magnetic field or an electric quadrupole in- 
teraction, then the neutron spin-flip produces a change 
of the ground state energy AE = Em — Em±i- This en- 
ergy change is transferred to the scattered neutron. The 
double differential scattering cross section^ is given by 
the following expressions: 

(1) 
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where a and a' are coherent and spin-incoherent scatter- 
ing lengths, W{Q) is the Debye- Waller factor and Eq is 
the incident neutron energy, S is the Dirac delta function. 
The long bar signs on top of equation ([T]) and ([2| mean 
averages of the quantity on top of which they stay. If 
the sample contains one type of isotope then — is 

zero. Also W 1 ± 



^ w I because AE is usually much 

less than the incident neutron energy Eq. In this case 
2/3 of incoherent scattering will be spin- flip scattering. 
The hyperfine splitting lies typically in the energy range 
of a few iJ^eV. The inelastic spin-flip scattering of neu- 
trons from the nuclear spins can yield this information 
provided the neutron spectrometer has the required res- 
olution of about IfieV or less and also the incoherent 
scattering of the nucleus is strong enough. Because of 
relatively large spin incoherent cross section of natural 
Nd, the Nd-based compounds are very much suitable for 
the studies of nuclear spin excitations. The Nd has the 
natural abundances of 12.18% and 8.29% of ^-^^Nd and 
^*^Nd isotopes, respectively. Both of these isotopes have 
nuclear spin of / = 7/2 and their incoherent scattering 
cross sectiondP^ are relatively large, 55 ±7 and 5±5 barn 
for ^''^Nd and ^''^Nd , respectively. Taking into account 
the natural abundance of these isotopes the total spin 
incoherent scattering of the natural Nd is CTi(spin)= 7.f 
barn. 

A Polycrystalline NdFeAsO sample was prepared as 
described befor^Sl, xhe lattice constants (a = 3.966, c 
= 8.589 A, V = 135.1 A3) were extracted from Rietveld 
fits against laboratory powder X-ray diffraction data and 
are in excellent agreement with our previous worU^^'^ 
and recently reported data on crushed single crystalJ^^. 
These fits also revealed minor Nd203 and FeAs impuri- 
ties with a total weight fraction less than 5%, which is 
comparable to Ref 15, and which evidently do not affect 
the lattice dimensions of the main phase. No evidence 
for the presence of Fe was found. About 5 g of pow- 
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FIG. 1: Energy spectra of NdFeAsO at several temperatures. 



der NdFeAsO sample thus prepared was placed inside a 
fiat Al sample holder and was attached to the sample 
stick of a standard '^He cryostat of the high-resolution 
back-scattering neutron spectrometer IN16 of the Insti- 
tut Laue-Langevin in Grenoble. The neutron wavelength 
was A = 6.27 A. The inelastic signals from the sample 
were measured at several temperatures in the tempera- 
ture range 1.7 — 25 K. In order to measure the inelas- 
tic signal from NdFeAsO at lower temperature down to 
about 55 mK the sample was then placed inside the annu- 
lar space of a Cu double- walled cylindrical sample holder 
and was attached to the '^He-'^IIe dilution insert of the 
cryostat. 

Fig. [T] shows the inelastic spectra from NdFeAsO sam- 
ple at several temperatures. The data for all Q has been 
integrated. We however checked carefully the Q depen- 
dence of the signals and convinced ourselves hat no such 
dependence exists. The measured Q range was from 
about 0.2 to 1.9 A~^. The lack of appreciable Q depen- 
dence of the inelastic signals due hyperfine interaction is 
expected and is checked by us in a variety of Nd-based 
compounds. At T = 0.055 K clear inelastic signals were 
seen at .E = 1.600 ± 0.003 /xeV on both energy loss and 



3 




FIG. 2: Temperature variation of the energy of the inelastic 
peak of NdFeAsO. 



energy gain sides. The inelastic signals move to the cen- 
tral elastic peak at higher temperature and finally almost 
merg with it at T/vd « 6 K. However the inelastic signals 
at a very small energy seem to persist at higher tempera- 
ture but could not be evaluated due to their proximity to 
the central elastic peak. At T = 25 K the signal is purely 
elastic and its width is explained by the resolution func- 
tion of the spectrometer. The spectra have been fitted 
to the convoluted resolution function of the spectrome- 
ter measured by the signal from the NdFeAsO sample at 
T = 25 K. Fig. [2] shows the temperature dependence of 
the energy of nuclear excitations. The fitting procedure 
does not work any more at T = 6 K due to the proxim- 
ity of the signal to the central elastic peak as we already 
noted before. The temperature variation of the energy of 
nuclear spin excitations shows that the Nd ordering pro- 
cess consists of two steps. At temperatures above 2 K the 
Nd ordering is perhaps induced by the increasing Nd-Fe 
exchange interaction whereas below 2 K the ordering is 
due to the onset of the Nd-Nd interaction. Unfortunately 
due to the limited neutron beam time we could not mea- 
sure the spectra at smaller temperature intervals which 



FIG. 3: Plot of ener gy of the inelastic signals in Nd2CuOj^, 
Nd metaP, NdCuJ^S, NdGaOM NdFeOj^, NdMgJ^ 
NdCoJ^, NdMnOiSl and NdAlJ^Il vs. the corresponding 
electronic magnetic moment of Nd in these compounds de- 
termined by the refinement of the magnetic structure using 
magnetic neutron diffraction intensities. We used the energy 
of the nuclear spin excitation determined in the present ex- 
periment to determine the magnetic moment of Nd ion in 
NdFeAsO at T = 0.055 K by using the hnear plot. The 
red filled square in the plot shows the corresponding data for 
NdFeAsO. The magnetic moment determined from the linear 
plot is 1.26/ifl. 



would establish the two stages of the Nd ordering pro- 
cess more clearly. At millikelvin temperature the effect 
of detailed balance is observed in the inelastic signals in 
the energy loss and energy gain sides. If one compares 
carefully the spectra at T = 0.055 K and T = 0.25 K 
in Fig. [T] then one will notice that at T = 0.055 K the 
inelastic signal in the energy loss side is weaker than that 
in the energy gain side whereas at T = 0.25 K the signals 
on both sides are about equal. 

We have studied hyperfine interactions in several Nd- 
based compound recently by the high resolution neutron 
scattering technique and have found a linear relationship 
between the energy of nuclear spin excitations at low tem- 
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perature and the ordered magnetic moment of the Nd 
ions in these compounds. One expects at low temper- 
ature in insulating Nd compounds an ordered magnetic 
moment of 3/is but due to the crystal-field effects one 
obtains a much lower ordered magnetic moments in neu- 
tron diffraction experiments. In Figure [3] w e show a plot 
of energy of inelastic peaks in Nd2CuOii^, Nd metaP^, 
NdCu22\, NdGaOa^^, NdFeOa^s, NdMg;^, NdCo:j23l 
NdMnOsj^and NdAlj^vs. the corresponding electronic 
magnetic moment of Nd ions in these compounds deter- 
mined by the refinement of the magnetic structure using 
magnetic neutron diffraction intensities. The data lie ap- 
proximately on a straight line showing that the hyperfine 
field at the nucleus is approximately proportional to the 
electronic magnetic moment. The slope of the linear fit 
of the data gives a value of 1.27zt0.03fj,eV / ^b- It is to be 
noted that the data for the hyperfine splitting is rather 
accurate whereas the magnetic moments determined by 
neutron diffraction have large standard deviations and 
are dependent on the magnetic structure models. The 
magnetic structures are seldom determined unambigu- 
ously and the magnetic moment determined from the re- 
finement of a magnetic structure model is relatively un- 
certain. In such cases the investigation of the low energy 
excitations described here can be of additional heljpS. 
This is specially true for the complex magnetic structures 
with two magnetic sub-lattices of which one sublattice 
contains Nd. Such complex magnetic structures, such as 



the parent compounds of newly discovered Fe-based su- 
perconductors, colossal magnetoresistive manganites and 
some multiferroic materials, are currently under intense 
study. Here we have used the nuclear spin excitation 
energy of NdFeAsO determined from the present exper- 
iment to determine the magnetic moment of Nd ion in 
NdFeAsO at T = 0.055 K by using the hnear plot shown 
in Fig. [3] The magnetic moment determined from the 
linear plot is 1.26/iB. Thus the magnetic moment of Nd 
ions are reduced considerably by the crystal field effect. 
The ordered magnetic moment of Nd in NdFeAsO has 
not yet been unambigusly reported by neutron diffrac- 
tion. The present estimation of the ordered magnetic 
moment of Nd in NdFeAsO remains to be checked by 
neutron diffraction experiment at very low temperature. 
Detailed neutron diffraction investigation on NdFeAsO 
single crystals is urgently needed. It would also be useful 
to study the crystal field effects in NdFeAsO by inelastic 
neutron scattering. 

In conclusion we have observed inelastic neutron scat- 
tering signals at low temperature in NdFeAsO, which 
give direct evidence for the magnetic ordering of Nd ions 
and we have estimated from the linear relationship es- 
tablished for Nd-based compounds the ordered magnetic 
moment of Nd to be 1.26 fiB- This estimate of the Nd 
magnetic moment is independent of the magnetic struc- 
ture model. 



^ Y. Kamikahara et al. J. Am. Chem. Soc. 130, 3296 (2008). 

^ D. Johrendt and R. Poettgen, Angew. Chem. Int. Ed 
47,4782 (2008). 

^ Z.-A. Ren et al., Chin. Phys. Lett. 25, 2215 (2008). 

" J.G. Bednorz and K.A. Mueller, Z. Phys. B: Condens. Mat- 
ter 64, 189 (1986). 

^ V. Johnson and W. Jeitschko, J. Solid State Chem. 11, 161 
(1974). 

^ J. Dong et al., Europhys. Lett. 83, 27006 (2008). 

G.F. Chen et al., Phys. Rev. Lett. 100, 247002 (2008). 
® T. Nomura et al., Supercond. Sci. Technol. 21, 125028 

(2008). 

^ C. de la Cruz et al.. Nature 453, 899 (2008). 

° W. Lu et al.. Solid State Comm. 148, 168 (2008). 

^ H.-H. Wen et al., Europhys. Lett. 82, 17009 (2008). 

^ P.J. Baker, S.R. Giblin, F.L. Pratt, R.H. Liu, G. Wu, 
X.H. Chen, M.J. Pitcher, D.R. Parker, S.J. Clarke and 
S.J. Blundell, New J. Phys. 11, 025010 (2009). 

^ Y. Qiu et al., Phys. Rev. Lett. 101, 257002 (2008). 

-'^ Y. Chen et al., Phys. Rev. B 78, 064515 (2008). 

W. Tian, W. Ratcliff II, M.G. Kim, J.-Q. Yan, P.A. Kien- 
zle, Q. Huang, B. Jensen, K.W. Dennis, R.W. McCal- 
lum, T.A. Lograsso, R.J. McQueeney, A.I. Goldman, J.W. 
Lynn, and A. Kreyssig, Phys. Rev. B 82, R060514 (2010). 

^ A. Marcinkova, D.A.M. Grist, ... S.J. Makowski, J. A. 
Rodgers, P.F. Henry, J. P. Attfield and J.W.G. Bos, Chem. 
Mater. 21, 264 (2009). 

A. Marcinkova, D.A.M. Grist, I. Margiolaki, T.C. Hansen, 



S. Margadonna and J.-W. G. Bos, Phys. Rev. B 81, 064511 

(2010) . 

J.-Q. Yan, Q. Xing, B. Jensen, H. Xu, K.W. Dennis, 
R.W. McCuUum, T.A. Lograsso, Phys. Rev. B 84, 012501 

(2011) . 

^® T. Chatterji and B. Frick, Physica B 276-278, 252 (2000) 
2° T. Chatterji and B. Frick, Appl. Phys. A 74[Suppl.], S652 
(2002) 

2^ T. Chatterji and B. Frick, Physica B 350, elll (2004) 
T. Chatterji and B. Frick, Sohd State Comm. 131, 453 
(2004). 

T. Chatterji, G.J. Schneider and R.M. Galera, Phys. Rev. 
B 78, 012411 (2008). 

T. Chatterji, G.J. Schneider, L. van Eijk, B. Frick and 
D. Bhattacharya, J. Phys: Condens. Matter 21, 126003 
(2009). 

T. Chatterji, G.J. Schneider and J. Persson, Phys. Rev. B 
79, 132408 (2009). 

A. Heidemann, Z. Phys. 238, 208 (1970) 
A. Heidemann and B. Alefeld, Neutron Inelastic Scatter- 
ing, International Atomic Energy agency, Vienna, 1972 
V.F. Sears in International Tables for Crystallography, Sec- 
ond Edition, vol. C, p. 445, ed. A. J.G. Wilson and E. 
Prince, Kluwer Academic Publishers (1999) 
^® R. Przenioslo, I. Sosnowska and B. Frick, J. Magn. Magn. 
Mat. 305, 186 (2006). 



